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Abstract 

On 2004 December 27, a giant flare from the soft gamma repeater 1806—20 was 
observed. The radiation mechanism of the initial peak of the flare would be con¬ 
troversial. In this letter we point out that very high-energy cosmic rays would be 
produced in the case that the flare was caused by internal shocks, as is usually con¬ 
sidered for gamma-ray bursts. The highest energy of cosmic rays can reach 10^® eV, 
if the Lorentz factor of the shocks is sufficiently high. Future observations of cosmic 
rays will inform us about the mechanism of the giant flare. 
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1. Introduction 

Soft gamma repeaters (SGRs) are X-/gamma-ray transient sources that show periods 
of bursting activity separated by long intervals of quiescence. They are galactic and LMG 
populations, and are considered to originate from neutron stars with intense (^ 10^® G) magnetic 
fields (magnetars) (Thompson, Duncan 2001). 

On 2004 December 27, a giant flare from SGR 1806—20 was observed by several detectors 
(Hurley et al. 2005; Terasawa et al. 2005; Palmer et ah 2005; Mazets et ah 2005). The initial 
peak of the giant flare had ~ 0.6 s duration. The isotropic-equivalent energy would have 
been E.^ ~ erg and the peak luminosity in the first 125 ms would have been L.y ~ 

lO'^^dis^ erg s“^, where dis = d/15 kpc, and d is an uncertain distance to the source (Gorbel, 
Eikenberry 2004; Gameron et al. 2005; McGlure-Griffiths, Gaensler 2005). 

In the context of the magnetar model (Thompson, Duncan 2001), the giant flare arose 
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from a hot expanding fireball, which was only weakly polluted by baryons. The initial optical 
depth for pair creation is extremely large in this model. As a result, this radiation-pair plasma 
expands relativistically. The Lorentz factor, 7 , increases in proportion to the radius, i?, and 
the temperature in the comoving frame decreases as T oc R~^ (Piran et ah 1993). The photons 
decouple when the temperature decreases below ~ 20 keV (Nakar et ah 2005). At this stage 
the Lorentz factor becomes ~ 10, and the radius becomes R ~ lOi^o, where Ro ~ 10 km is 
the size of the magnetar. The photon temperature that we observe is the initial temperature 
independently of 7 at the decouple stage because of the relativistic blue-shift. The predicted 
photon spectrum from this “pure radiation-pair hreball model” is, therefore, a quasi-blackbody 
with a temperature of T ~ 200 keV (Hurley et ah 2005). The time-resolved (125 ms) energy 
spectrum from the RHESSI particle detector is consistent with that of a black body whose 
temperature is ~ 200 keV, though the main RHESSI spectroscopy detectors were saturated 
during the peak (Boggs et al. 2005). 

However, there are so many ambiguities that one cannot conhrm the above picture. The 
spectrum of the initial peak of the giant flare is not yet well determined. Hurley et al. (2005) 
reported on the cooling blackbody spectrum, which is consistent with the pure radiation-pair 
hreball model. However, Mazets et al. (2005) derived a power-law spectral shape (a = —0.7) 
with an exponential cut-off at 800 keV. There was another independent measurement, resulting 
in a power-law spectrum {a = —0.2) with an exponential cut-off at 480 keV (Palmer et al. 2005). 
Furthermore, another giant hare in SGR 0526—66 may also have had a nonthermal spectrum 
(Fenimore et al. 1996). Therefore, we still have no dehnite evidence to interpret that the hare 
on 2004 December 27 was a thermal one. 

Another ambiguous point is the initial speed of the outhow or, almost equivalently, the 
baryon richness. The radio afterglow detected after the giant hare may give us some hints to 
solve this problem (Cameron et al. 2005; Gaensler et al. 2005; Taylor et al. 2005). Wang et 
al. (2005) showed that the same mechanism as established for gamma-ray burst (GRB) after¬ 
glows (decelerating-outhow model) can explain the radio afterglow of the giant hare, suggesting 
that the initial outhow was highly relativistic. On the other hand, Granot et al. (2005) and 
Gelfand et al. (2005) succeeded to ht the radio afterglow by their model, in which the outhow 
is initially non-relativistic and contains baryonic material of more mass than the 10 ^^ g ejected 
from the magnetar, itself. The baryon amount in this picture, at hrst glance, may contra¬ 
dict the optically thin emission of the hare. Therefore, some authors (Eichler 2005; Dai et 
al. 2005) propose multi-component models, in which relativistic (baryon-poor) components and 
non-relativistic (baryon-rich) components outhow from the surface. 

Jet collimation has been suggested by Yamazaki et al. (2005) (Y05) through reproducing 
the observed light curve of the initial peak observed by GEOTAIL with a 5.48 ms time resolu¬ 
tion (Terasawa et al. 2005). Using a simple emission model for relativistically moving matter 
(Yamazaki et al. 2003), they derived an upper limit of the jet opening half-angle of 0.2 rad. In 
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this model the initial peak in the giant flare arises from internal shocks in relativistic jets. The 
radius, i?i, where the shock front starts to emit photons, is estimated to be Ri = 2.6 x 10 ® 7 ^ cm, 
which is much larger than the radius, 7 i?o ~ 10 ^ cm, derived from the pure radiation-pair fireball 
model. Considering the opening half-angle of the jet, < 8/7 (Y05), the collimation-corrected 
energy of the jet may be 10^® erg if the radiation efficiency is assumed to be 10%. The observed 
proper motion of the centroid of a radio image (Taylor et ah 2005) may support the collimated 
outflow. The relativistic jet is significantly decelerated due to the sideways expansion within 
ten minutes after the giant flare (Rhoads 1999; Sari et ah 1999; Yamazaki et ah 2005), so that 
the radio image expands non-relativistically (Taylor et ah 2005). Introducing a non-relativistic 
component besides the relativistic jet, the observed expansion law of the radio image may also 
be explained (Dai et ah 2005). 

Therefore, there are three controversial questions concerning this event: 1) Are there 
baryonic flows during the initial stage? 2 ) Are there relativistic components in the baryonic 
flows? 3) Does the flow have a jet-collimated structure? We cannot, at present, definitely 
answer these questions. In this paper we show that observations of ultra-high-energy cosmic 
rays (UHECRs) or neutrinos may resolve the above-mentioned problems (Eichler 2003; loka 
et al. 2005). loka et al. (2005) extensively discussed the possibility of photopion production 
and the resultant neutrino detection from this event. On the other hand, our discussion places 
emphasis on the detection possibility of UHECRs delayed from the giant flare, based on the 
model in Y05, which can explain the light curve of this event. 

If there is a shocked outflow of baryons, the shock waves may accelerate non-thermal 
protons to very high energies. As for extragalactic GRBs, Waxman (1995); Vietri (1995) 
pointed out that GRB is a candidate of UHECR sources, and Milgrom, Usov (1995) discussed 
possible observations of UHECRs. In this paper we discuss the possibility of UHECR production 
in the giant flare of SGR 1806—20. Combining neutrino observations, UHECR detections may 
provide information on the above questions 1) and 2). In section 2 we obtain the maximum 
energy of UHECRs accelerated in shock waves based on the jet model of Y05. In section 3 the 
expected number flux of UHECRs, which depends on the structure of the galactic magnetic 
fields, is discussed. Our conclusions are summarized in section 4. 

2. Maximum Energy of Cosmic Rays 

In this section we adopt the jet model of Y05 and estimate the maximum energy of 
cosmic rays that we can observe. We expect that even another internal shock model may 
predict similar maximum energies to ours, as long as the assumed energy and time-scale of the 
flare are close to ours (see e.g. Eichler 2005). 

In the model of Y05 a shell moving with the Lorentz factor, 7 , emits gamma rays from 
R = R[ = 2.6 X 10 ® 7 ^ cm to Re — 10i?i. The duration in the comoving frame is {Re — Ri) /( 07 ) ~ 
Re/{c^), which corresponds to the shell crossing time scale. Therefore, the shell width in the 
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observer frame is about i?e/ 7 ^, while the shell width used conventionally in the GRB models is 
Rill'll which implies Rc ~ Ri. Thus, the shell width obtained from the GEOTAIL light curve 
is relatively thick. As a result, the physical conditions around R^ R\ and R^ Re are different. 

It is difficult to estimate the magnetic field in the shell, because the light curve pro¬ 
vides information only on the kinematic properties. The plasma in the shell coming from the 
magnetar may be strongly magnetized. A part of the energy of the magnetic field may be 
expended in the acceleration phase, while various plasma instabilities may enhance the held, as 
is discussed by Medvedev, Loeb (1999) for GRBs. Though we recognize the above difficulties, 
we assume that the energy density of the magnetic held is proportional to the photon energy 
density. The luminosity, which is Lorentz invariant, is E^c/Re = 1.2 x IO'^^E' 467 ^^ erg s“^, where 
= 10 ^®E 46 erg. In this case the photon energy density in the comoving frame is approximated 
as U^ = E^/^AnR^Re). From R^/(87r) = £ 3 ^ 7 , we obtain R(i?) = 1.1 x 

G, where 7 = IOO 72 . The energy ratio, £b, can be larger than unity, if the baryonic energy 
density is much larger than which implies that the gamma-ray emission is an inefficient 
process. 

The shocks may accelerate non-thermal protons. The maximum energy of accelerated 
particles, Emax, is restricted by the condition that the particle Larmor radius, = E/jeB, is 
smaller than the size scale of the emitting region, {R — Ri)/'y ~ R/l- This condition gives the 
maximum energy, El, max = 8.3 x eV, which is independent of R. 

Another condition to generate cosmic rays is that the cooling time scale of protons should 
be longer than the dynamical time scale Tdyn = R/{c'y). As long as we consider energies below 
El, max, the acceleration time scale (~rL/c) may be shorter than Tdyn. If the protons accelerated 
to high energies cool down before they escape from the shell, those particles cannot become 
cosmic rays. For proton-synchrotron cooling we obtain the maximum energy of E^^m max = 
4.6 X 10 ‘^^ e ^^- flE ^ i { R/RO eV. 

We should comment on the “adiabatic cooling” due to shell expansion. While the shell 
coasts from E, to Re, the shocked region grows via shock propagation. The adiabatic invariance, 
Br^, may lead to the cooling of particles. However, we are not snre whether the adiabatic 
condition is satished or not for a disturbed magnetic held in the shocked region. In any case, 
we assume that particles escape from the shell on the dynamical timescale at each radius via 
stochastic processes. Therefore, we neglect adiabatic cooling in our case. 

The maximum energy of cosmic rays is determined by min(EL,max,Esyn,max)- The results 
are plotted in hgure 1. For a smaller R the maximum energy is determined by Esyn max , because 
the magnetic held becomes stronger. As R increases the maximum energy is determined by 
EL,max, rather than E^yn max , because the weaker magnetic held cannot conhne high-energy 
particles in the shell. For 7 ^ 30 the radius where the shock occurs (oc 7 ^) is smaller, so that 
the energy limit is determined by Esyn,max only between R^ and Eg- On the other hand, for 
7 ^ 40, the limit is determined by only El, max- The values, Emax, achieved in shocked shells 
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Fig. 1. Maximum energy, Fmax, determined by i?L,max and £^syn,max for £6 = 1 and = 1. The Lorentz 
factor, 7 , is assumed to be (a) 10, (b) 25, (c) 30, (d) 40, (e) 100, and (f) 200. 

with 7 are the maximum values in each line in hgure 1. For £b = 1 and = 1, i^max niay be 

a few times 10^® eV. 

Another cooling process that we should take into account is photopion creation. If the 
proton photopion “optical depth” is much higher than unity, protons lose their energies in the 
flare source before they escape from the emitting region (Asano, Takahara 2003; Asano 2005). 
The condition to create pions is Ec^e^y > 0.2 GeV^, where F^cr and €j are the energies of a 
nucleon and an interacting photon, respectively. In order to estimate the energy loss rate due 
to photopion creation, we need the soft photon spectrum in the flare, which is unknown. 

If the giant flare is thermal emission of ~ 200 keV, the number of soft photons that 
interact with protons of ~ 10^® eV is too small to create pions on the dynamical time scale. 
However, it may be premature to conclude that the spectrum observed by RHESSI is thermal. 
As suggested in Mazets et ah (2005) and Palmer et ah (2005), there is a possibility that soft 
photons due to non-thermal electrons may distribute below ~ 200 keV with the spectrum 
n{ery) oc e“^, as observed in standard GRBs. 

As the most pessimistic case, we assume that power-law photons with p = 1 or 1.5 
dominate below 200 keV without the low-energy cut-off. Above 200 keV we adopt the Planck 
spectrum with 200 keV, though the spectrum shape above 200 keV is not important for the 
cooling process of ultra-high-energy particles. Although the number of power-law photons is 
much larger than the high-energy photons, the energy contribution of the soft photons is not 
important for p <2. Using the same cross section and method as Asano (2005), we estimate 
the time scale of the photopion creation, Tp, near i?max- The results are plotted in hgure 2. 

The low-energy bump appearing in hgure 2. is due to the “thermal” photons above 200 
keV. Above this energy nucleons interact with the power-law photons. In this case the number 
of interacting photons {n{ery)dery, ~ 0 . 27 ^GeV^/E'cR) is independent of Ecr for p = 1. As a 
result, the time scale is nearly constant in this energy region. For p = 1.5 the time scale becomes 
shorter as Eqr increases (Tp, oc If the power-law spectrum has a low-energy cut-oh due 
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Fig. 2. Time scale ratio Tdyn/Tpi vs. energy of nucleons for photon index p = 1 (solid) and p = 1.5 
(dashed). The model parameters ( 7 ,i?/i?i) are (a) (10,10), (b) (25,10), (c) (40,1), and (d) (100,1). The 
total photon energy, Ej, is fixed as 10"^® erg. The dashed lines of (b) and (c) nearly overlap. 

to synchrotron self-absorption etc., the time scale will increase from the corresponding energy. 

For 7 = 10 the photopion prodnction is crucial, and resultant high-energy neutrinos 
will be emitted as is the case in loka et ah (2005). However, if 7 is large enough to generate 
UHECRs, the photopion production is inefficient. Even for 7 = 25, UHECRs will survive and 
neutrinos may not be observed, unless the photon spectra are extremely soft, such as p = 1.5. 

Unless the energy of soft photons below 50 keV is much larger than the observed energy 
by RHESSI, there is a possibility that UHECRs of > 10^® eV would come from the SGR 
1806—20 giant flare. At present, considering the ambiguity in the parameter values, we cannot 
completely preclude the detection possibility of UHECRs above 10^° eV (for example, Eb = 10, 
7 = 30, and = 3). 

3. Propagation of UHECRs 

The source location of the flare is about 10° from the galactic center. In the region 
around the galactic center the magnetic field is highly uncertain. If the magnetic fields are 
well represented by regular helds along the spiral arms, the time delay may be on the order of 
~ 10 yrs, even for 10^° eV (Alvarez-Muniz et ah 2002). Turbulent magnetic helds, whose scale 
is 10-100 pc (Beck 2001), may shorten the above estimate. In the most pessimistic case, the 
time delay becomes on the order of thousands of yrs (Eichler 2005). 

If the energy of cosmic rays above 10^® eV is Pig % of the total photon energy ~ 10^® 
erg, 2PiQdi ^particles are detected per one year by 1000 km^ detectors, such as AUGER 
or Telescope Array, where AT 3 is the dispersion of the timedelay normalized by 1000 yr. If the 
particles are conhned to 0.1 rad of the galactic plane, the above detection rate can be sufficient 
signals, even for AT 3 = 1 (Eichler 2005). 
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4. Conclusions and Discussion 


If the giant flare from SGR 1806—20 on 2004 December 27 is due to internal shocks 
in relativistic jets, we found that the shocks may produce UHECRs of up to rsj 10^9 eV. The 
maximum energy is similar to an estimate by Eichler (2005). In standard GRBs, on the other 
hand, UHEGRs are unlikely to be observed, since they lose energy before escaping from the 
shell via photopion production (Asano 2005). Even if UHEGRs are produced in GRBs, the 
time delay between cosmological bursts and UHEGRs is too large to establish any connection 
between the GRBs and UHEGRs. On the other hand, because the giant flare of SGRs is less 
luminous than standard GRBs, UHEGRs can escape from the shell without losing their energy. 

In our baryon-loaded jet model, if the bulk Lorentz factor of the outflow is high enough, 
neutrinos that are produced from the decay of charged pions may not come from this flare. 
This is consistent with the results of loka et ah (2005). On the other hand, if the outflow 
is non-relativistic, neutrinos could be produced via p-p collisions, which is neglected in our 
calculation, in addition to photopion production. Hence, UHEGRs and high-energy neutrino 
observations become diagnostic tools to investigate the properties of the outflows at an early 
phase. There are three possible cases: A) detection of neutrinos, B) detection of UHEGRs, 
but no neutrinos, and G) no detection of high-energy particles. We may reach the following 
conclusion irrelevant of the models. Gases A and B mean there are baryons in the flare stage, 
but the Lorentz factor is not very large in case A (see also Gelfand et ah 2005). Gase G 
implies that there are negligible baryons (pure radiation-pair fireball model), or the UHEGR 
production efficiency is low. 

Because the arrival time of UHEGRs depends on a number of factors, it cannot be 
predicted exactly. It might be at present that we will detect UHEGRs coming from the other 
two SGRs (SGR 0526—66, SGR 1900—14) that previously caused giant flares. Giant flares of 
SGRs could produce a large amount of UHEGRs, which may explain the origin of doublet and/or 
triplet events, though there is no evidence of any correlation between such events observed by 
AGASA and SGRs so far. Anomalous X-ray pulsars (AXPs) may be the same kind of objects as 
SGRs and older, but less active than SGRs (Kulkarni et al. 2003). They might have produced 
UHEGRs via giant flares thousands of years ago. We could detect signals from past activities 
of AXPs by UHEGR observations. 


This work was supported in part by JSPS Research Fellowship for Young Scientists (RY). 
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